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Abstract. A comparison between the two tracers of magnetic field mirror asym- 
metry in solar active regions, twist and current helicity, is presented. It is shown 
that for individual active regions these tracers do not possess visible similarity while 
averaging by time over the solar cycle, or by latitude, reveals similarities in their 
behaviour. The main property of the dataset is anti-symmetry over the solar equator. 
Considering the evolution of helical properties over the solar cycle we find signatures 
of a possible sign change at the beginning of the cycle, though more systematic 
observational data are required for a definite confirmation. We discuss the role of 
both tracers in the context of the solar dynamo theory. 
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1. Introduction 



Contemporaty astronomical observations suggest two proxies for mirror 
asymmetry in solar active regions, namely, the current helicity and 
twist of magnetic field. Both proxies, averaged over a suitable part of 
active regions, have been measured recently for several hundred active 
regions over more than one solar cycle (see, e.g., Bao and Zhang, 1998 
and references therein). Both proxies demonstrate, to some extent, an 
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antisymmetry with respect to the solar equator as well as a cyclic 
behaviour on timescales of the solar activity cycle as traced by sunspots. 

Solar magnetic field structure is complicated enough to allow many 
proxies for its mirror asymmetry, which may not necessarily be pro- 
portional one to another. As is natural to expect, current helicity and 
twist have no simple relation between each other. 

A detailed comparison between current helicity and twist as two 
proxies of the mirror asymmetry of the solar magnetic field is interest- 
ing in the context of solar dynamo theory; the key driver of the solar 
dynamo suggested by Parker (1955) is the a-cffcct originating in the 
mirror asymmetry of solar convection and magnetic fields. Proxies of 
mirror asymmetry of the magnetic field in solar active regions provide 
a unique observational approach for the direct verification and obser- 
vation of the a-cffcct. The link between the a-effcct and the proxies 
under discussion is usually given in terms of the current helicity {e.g., 
Kleeorin et al., 2003). It would however be risky to insist that solar 
dynamo theory is well-developed enough to disregard twist as an alter- 
native tracer of the a-cffect. A conventional mean-field description of 
the a-effect deals with quantities averaged over substantial temporal 
or spatial domains rather than over an individual active region. If the 
similarity of current helicity and twist as tracers of mirror asymmetry 
becomes more pronounced after averaging, it means that both proxies 
are reasonable tracers of the a-effect and support the conventional 
concept of the solar dynamo. 

In the present paper we compare current helicity and twist data. 
First, we study the proxies for an individual active region to demon- 
strate that the correlation between current helicity and twist is rather 
weak. Then we demonstrate that the similarity becomes much more 
pronounced after temporal or spatial averaging. 

Our investigation is based on data obtained at the Huairou Solar 
Observing station of the National Observatories of China (Zhang and 
Bao, 1998). A previous analysis of current helicity and twist based on 
this data was presented by Zhang et al. (2002), see also Kuzanyan et 
al. (2000). Here we use a larger data set and improve the statistical 
analysis, as well as embedding the study in the context of the solar 
dynamo in a more explicit form. 

The paper is organized as follows. We briefly review the concepts of 
current helicity and twist as they are exploited in theoretical studies 
and give their observational proxies (Section 2). Then we describe the 
observational data set (Section 3), compare the current helicity and 
twist for a particular active region (Section 4) and then present the 
data after temporal or spatial averaging (Section 5). Section 6 contains 
a more detailed analysis of the antisymmetry of the proxies with respect 
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to the solar equator. We discuss the consequences for solar dynamo 
theory from this analysis in Section 7. 



2. Current Helicity and Twist and the a-effect 

The conventional parametrization of the magnetic contribution to the 
a-effect (e.g. Kleeorin and Rogachevskii, 1999) is based on the current 
helicity x'^ =<t» • j>, where j = curlb is the electric current and b is 
the (small-scale) magnetic field. Because curlb is calculated from the 
surface magnetic field distribution, the only electric current component 
that can be derived is (curlb);^. As a consequence of these restrictions, 
the observable quantity is 

H, = {b,icuAb),) , (1) 

where x,y,z are local cartesian coordinates connected with a point on 
the solar surface, and the z-axis is normal to the surface (Bao and 
Zhang, 1998; see also Abramenko et al, 1997; Pevtsov et al, 1994). In 
the framework of the hypothesis of local homogeneity and isotropy this 
value is 1/3 of the current helicity x^- 

Observations {e.g. Zhang, and Bao, 1998) provide another proxy for 
the mirror asymmetry of the magnetic field, i.e. twist (Woltjer, 1958) 
which comes from studies of magnetic fields in the solar atmosphere, 
where conductivity is high. However, because of low-beta condition the 
magnetic field can be described as force-free. Furthermore, according to 
Maxwell's equations, the magnetic field is a Beltrami field, i.e., curlb = 
affb, where the parameter as is the twist. In the solar interior, however, 
the magnetic field is not a Beltrami field and twist can be understood 
as as =< curlb ■ b/6^>. Of course, this definition does not coincide 
with that of the current helicity. The observational equivalent of the 
quantity as is the ratio <jz/bz>- The notation as is generally used for 
twist in the solar physics literature, though it seems rather confusing 
from the viewpoint of solar dynamo theory. The details of calculation 
of twist and helicity from magnetographic observational data are given 
in literature (e.g. Wang et al., 1996; Bao and Zhang, 1998; Zhang and 
Bao, 1998). 



3. Observational Data 

The observational data used in our analysis were obtained at the Huai- 
rou Solar Observing station of the National Astronomical Observatories 
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of China. A magnetograph using the Fel 5324 A spectral hne deter- 
mines the magnetic field values at the level of photosphere. The data 
are obtained using a CCD camera with 512 x 512 pixels over the whole 
magnetogram. The entire image size is comparable with the size of an 
active region, which at about 2 x 10^ m is comparable with the depth 
of the solar convective zone. 

However, the observational technique allows the line-of-sight field 
component to be determined with a much higher precision than 
the transverse components {bx and by). There are a number of other 
observational difficulties, such as resolving the 180° ambiguity in the di- 
rection of transverse field etc. The observational technique is described 
in detail by Wang et al. (1996), see also Abramenko et al. (1996). 

The observations arc restricted to active regions on the solar surface 
and we obtain information concerning the surface magnetic field and 
helicity only. Monitoring solar active regions while they are passing near 
the central meridian of the solar disc enables observers to determine 
the full surface magnetic field vector. The observed magnetic field is 
subjected to further analysis to determine the value of curlb. 

An observational programme to reveal the values of the twist and 
the current helicity density over the solar surface requires a systematic 
approach, both to the monitoring of magnetic fields in active regions 
and to the data reduction, in order to reduce the impact of noise. 
This work has been carried out by a number of research groups {e.g.. 
Seehafer 1990; Pevtsov et al. 1994; Rust and Kumar 1996; Abramenko 
et al. 1997; Bao and Zhang 1998; Kuzanyan et al. 2000). 

The present work analyzes the two systematic datasets of active 
regions. The first one consists of 422 active regions over the 10 years 
1988-1996 (Bao and Zhang, 1998). It has been used for theoretical 
analysis and further data reduction by Kuzanyan et al. (2000), Zhang 
et al. (2002), Kleeorin et al. (2003), Zhang et al. (2006) and Sokoloff et 
al. (2006). 

We also analyse a dataset which covers the three years at the begin- 
ning of the solar cycle 23, namely 1998-2000. This dataset was discussed 
earlier by Bao et al. (2000, 2002), and contains data for 64 active regions 
for which all the helical quantities were determined. The new data 
are obtained using the same technique and are processed in much the 
same way, as the earlier dataset of Bao and Zhang (1998) covering the 
ten year period 1988-1997, see also Zhang and Bao (1998). Thereafter, 
following Sokoloff et al. (2006), we merge these two sets of data and 
henceforth will consider them as a single continuous dataset of 486 
active regions. 

Observational work is ongoing, and much more data are due to be 
processed shortly {e.g., Xu et al., 2007). 
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Figure 1. Cumulative distribution function for the current helicity (upper panel) 
and twist (lower panel). The first coordinate of the point is the value of the pa- 
rameter under investigation for a particular active region. The second coordinate 
is the expected value of the standard deviations for a Gaussian quantity with the 
same mean and standard deviation which gives the same probability (see in text 
for details). For a Gaussian quantity, we must obtain straight lines, which are also 
shown on both panels. 



4. Helicity and Twist for a Particular Active Region 



We consider both observable parameters for a particular active region as 
random quantities. Ttie standard way to represent a random quantity is 
as a histogram presenting the percentage of data lying in a given range. 
We tried this method and found that the result is not very informative. 
The point is that our data set is small and the data are quite noisy. 
So we used a more advanced method and calculated a cummulative 
distribution function (c.d.f.) which is much more robust than the usual 
probability distribution function. 
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Figure 2. Scatter diagram for the twist and helicity normalized to the corresponding 
standard deviations. A correlation between the twist and helicity means that the 
distribution of points will form an ellipsoid, whose inclination to the coordinate axes 
provides the correlation coefficient. In the case under discussion, the correlation is 
seen to be weak, though a trend (shown by a dashed line) is noticeable. 

Let our set contain A'^ active regions. For the sake of definiteness 
consider a given value oi x = (or x'^ if appropriate). Let n active 
regions have twist lower than x. Then the probability for to be lower 
than X is estimated as P = n/N. Let ^ be a Gaussian variable with the 
same mean value /i and standard deviation a as ag and y the value 
for which the probability for — /i)/cj to be lower than y is P. The 
results for various x are shown by dots in coordinates (Figure 1; 

upper pannel is for twist and the lower one is for currrent helicity) 
and can be compared with the c.d.f. for a Gaussian distribution shown 
by the solid line. We see that the dots (twist) on the upper pannel 
are substantially closer to the straight line than on the lower panel 
(helicity). We conclude that twist is much closer to a Gaussian random 
quantity than helicity. Note that the link between both quantities is 
very nonlinear, and so at least one of the quantities (here the helicity) 
has to be non-Gaussian. 

Next, the twist and helicity of the observational data normalized 
to their means (/ig and fic, correspondingly) and standard deviations 
(cjg and fJc, correspondingly) are shown on a scatter diagram (Fig- 
ure 2). If the statistical dependence between the twist and helicity 
is strong, the cluster of points obtained must form an elongated el- 
lipse. In our case the axes of the ellipse are nearly parallel to the 
coordinate axes, and a weak correlation can be described by the re- 
lation (x'^ — fJ'c)/(^c = 0.006 + 0.1 (og — ^g)/(Tg obtained by the least 
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Figure 3. Helicity (upper panel) and twist (lower panel) for individual active regions 
as a function of the phase cycle (for the southern hemisphere). The cycle is not 
pronounced though one can see some predominance of negative values in both the 
helicity parameters, in accord with the polarity law. 



square method. The resulting correlation may be formally significant 
however this is difficult to confirm because of the non-Gaussian dis- 
tribution of helicity. Anyway, the correlation revealed appears to be 
robust with respect to discarding data that strongly deviates from the 
mean. Note, however, that theoretical considerations based on such 
weak correlations are not very reliable. 

For the helicity and twist data obtained over individual active re- 
gions, the solar cycle is not well-pronounced. As an example, in Figure 3 
we show helicity (upper panel) and twist (lower panel) as functions of 
the cycle phase for active regions in the southern hemisphere. The 
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periodic behaviour is hardly visible in both cases. Notice that separate 
presentation of the data for northern and southern hemisphere is re- 
quired because of the hemispheric rule according to which both tracers 
tend to have opposite signs in the northern and southern hemispheres. 

On the other hand, appropriate averaging of the data taken over 
rather narrow temporal or latitudinal intervals indicates similar cyclic 
behaviour for both quantities. This is described in the following section. 

5. Evolution of the Mean Values Over an Activity Cycle 

The mean values of the twist and helicity calculated over relatively 
narrow time or latitudinal intervals behave quite differently. Follow- 
ing Kleeorin et al. (2003), we divide the entire data set into two-year 
time intervals and plot the mean helicity (Figure 4, upper panel) and 
mean twist (Figure 4, lower panel) for each interval separately for 
the northern and southern hemisphere. The error boxes are calculated 
assuming the quantities are Gaussian. Figure 4 shows that the cyclic 
variations of the helicity and twist are very similar. Both parameters 
increase in absolute value in the middle of the cycle and decrease at 
the beginning and end of the cycle. Both helicity and twist change sign 
from one hemisphere to another. For both tracers, the cycle is seen 
more distinctly in the southern hemisphere; the cycle in the northern 
hemisphere is somewhat better pronounced in helicity than in twist. 
On the other hand, the mean values contain significant uncertainties. 
Therefore, we cannot confirm the hypothesis that the mirror asymmetry 
of the magnetic field changes sign in the course of an activity cycle as 
suggested by Hagino and Sakurai (2005). We have checked that the 
result is robust with respect to discarding data that strongly deviates 
from the mean. 

Note that, in comparison with Kleeorin et al. (2003), we use an 
extended database covering a longer period of observations. The obser- 
vational database and the procedure providing a synthetic description 
of a full cycle using partial observations of two successive cycles are 
given in Sokoloff et al. (2006). 

A similar picture was obtained by averaging the data over 5° - 
latitudinal bands (Figure 5). Again, both tracers display a mirror asym- 
metry with respect to the solar equator. The helicity data in the south- 
ern hemisphere are more regular than in the northern one, while the 
twist data, on the contrary, seem to be more regular in the northern 
hemisphere. Also, we do not see polarity inversions over latitudinal 
bands in either hemisphere. (As shown by the confidence intervals, 
the apparent polarity reversal at —30° for helicity is insignificant and 
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Figure 4- Two-year averaged mean helicity (upper panel) and twist (lower panel) 
values as a function of the cycle phase calculated separately for the northern and 
southern hemispheres. Blue squares correspond to the northern hemisphere while 
red diamonds correspond to the southern one. The data points averaged over both 
hemispheres, when plotting, have been artificially shifted in opposite sides for better 
performance. One can easily see the polarity law and evolution of both the values 
over the cycle. 



disappears if strongly deviating values are discarded). Notice that a 
similar result was obtained by Zhang et al. (2002) using a smaller 
dataset, see their Figure 4. 

In general, we conclude that the twist data averaged over time or 
latitude intervals can be used to determine quite reliably the behaviour 
of the helicity and vice versa. 

Isolating the particular time intervals and latitudinal bands and 
separating the data by hemispheres, we decrease significantly the num- 
ber of the values to be averaged in each case. Therefore, the available 
data appears to be insufficient for further fragmentation. In particular, 
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Figure 5. The helicity (upper panel) and twist (lower panel) values averaged over 
5° -latitudinal bands. One can easily see the polarity law. 



butterfly diagrams for the helicity (Sokoloff et al, 2006) and twist 
(Figure 6) based on these data are merely illustrative and cannot be 
used directly to draw conclusions on, say, the inversion of the sign of 
helicity. 

Note that a weak correlation between the helicity and twist data 
for individual active regions discussed in the previous section needs 
an explanation in the context of the revealed pronounced correlation 
between their mean values. We may suggest that there is significant 
correlation between the helicity and twist for individual active regions 
calculated in a particular latitudinal range and cycle phase, but this 
may be disguised by the general cyclic dependence of the entire dataset. 
Thouth the limited size of the available dataset does not allow us to 
verify this suggestion. 
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Figure 6. Butterfly diagrams for the helicity (upper panel) and twist (lower panel). 
The red crosses mark the positions of the active regions with negative values of the 
corresponding parameters, and the blue triangles correspond to the positive values. 
The two diagrams look very similar, and so can hardly be distinguished by sight. 



6. Statistics of Active Regions Violating the Polarity Law 

Of course, the polarity law for the hehcity and twist is not strictly 
fohowed. As shown by the helicity and twist measurements in individual 
active regions, there are a lot of active regions that violate this law. It 
turns out (Sokoloff et al., 2006) that the active regions in which the 
current helicity does not obey the polarity law are most frequently 
observed at the beginning of the cycle (cf. Tang and Le, 2005). Table 1 
presents the corresponding statistics both for the helicity and twist. 
The tendency of the law-breaking active regions to appear preferably 
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Table I. Statistics of active regions breaking 
the polarity law for the current helicity (n_ 
is the number of active regions before the 
cycle phase t and p is the relative number 
(probability) of the law-breaking active re- 
gions) and the twist {N- is the active region 
number and p is the probability). 



t 


n- 


N- 


P 


P 


0.18 


18 


16 


54 ± 8% 


48 ± 18% 


0.30 


60 


54 


46 ± 4% 


42 ± 9% 


0.43 


85 


91 


37 ± 3% 


40 ± 6% 


0.55 


101 


123 


32 ± 2% 


38 ± 5% 


0.68 


112 


144 


28 ± 2% 


36 ± 5% 


0.80 


121 


162 


26 ± 2% 


35 ± 4% 



at the beginning of the cycle is noticeable in the case of the twist as 
well, but this is much weaker than in the case of the helicity (in contrast 
to what might be expected according to Choudhuri et al. (2004). This 
raises a problem which challenges further development of the dynamo 
theory. 



7. Discussion 

In this paper we have shown that the behaviour of current helicity 
and twist is similar. Both these quantities are obtained from the same 
distribution of local values of vertical magnetic field hz and the electric 
current j^, properly averaged upon calculation. Therefore this similarity 
is expected, however not evident in advance. Although one can suggest 
other quantities to be calculated from the same vector magnetographic 
dataset, such as atest) which can also be used as a tracer of mirror 
asymmetry of magnetic fields (see, e.g. Hagino and Sakurai, 2005). We 
may suggest to compare these related quantities in the forthcoming 
papers. 

The above analysis exploits sunspots as tracers of processes in the 
region of magnetic field generation, i.e. in the region of dynamo action. 
Let us discuss the applicability of this approach. At the photospheric 
level magnetic pressure in sunspots is likely much larger than the gas 
pressure, and so the magnetic field can be considered in the vacuum 
approximation. This supports the application of potential magnetic 
field models for the solar corona and to some extent for the photosphere. 
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Correspondingly, we use twist as a quantity reflecting mirror asym- 
metry of magnetic field in the potential approximation. The situation 
becomes quite different in the sub-photospheric region as the gas pres- 
sure and kinetic energy become comparable with the magnetic energy 
just below the photosphere. Correspondingly, we use current helicity as 
a quantity reflecting the mirror asymmetry in the solar interior. 

The extent to which the current helicity data taken at the photos- 
pheric level represent values for the dynamo region seems to be much 
more delicate. The point is that time-distance helioseismology demon- 
strated that sunspots evanesces at the depth of 5-10 Mm only. One 
might conclude that this fact precludes using sunspots as tracers of 
physical processes in the solar interior. We believe that such an opin- 
ion woTild be an exaggeration. First of all, helioseismology gives the 
depth of the region with a temperature depression only. Below this 
depth a temperature excess is expected (Pononiarenko, 1972a,b; Parker, 
1974, 1976). It is reasonable to believe that the magnetic field below 
a sunspot cannot suppress convection even though it is not negligible. 
Note that according to modern understanding in the framework of a 
cluster model, the motion of magnetic flux tubes is rather independent 
of convective motion. This viewpoint is supported by direct observa- 
tions (Zhao et al., 2004; Gizon et al, 2003). In addition, long term 
monitoring of sunspot rotation demonstrates a clear solid-body com- 
ponent typical for tachocline motion (Ivanov, 2004). One may expect 
that sunspots mimic somehow the structure of the large-scale toroidal 
magnetic fields in the solar interior. Therefore, we may use the data on 
current helicity and twist as tracers of the dynamo mechanism, though 
this question requires further clarification. 
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